ARITHMETIC PROPERTIES OF PICARD-FUCHS EQUATIONS AND 

HOLONOMIC RECURRENCES 

ZANE KUN LI AND ALEXANDER W. WALKER 

f/-\ ' Abstract. The coefficient series of the holomorphic Picard-Fuchs differential equation associated 

with the periods of elliptic curves often have surprising number-theoretic properties. These have 
been widely studied in the case of the torsion- free, genus zero congruence subgroups of index 6 and 12 
(e.g. the Beauville families). Here, we consider arithmetic properties of the Picard-Fuchs solutions 
associated to general elliptic families, with a particular focus on the index 24 congruence subgroups. 
We prove that elliptic families with rational parameters admit linear reparametrizations such that 
their associated Picard-Fuchs solutions lie in Z[[t]]. A sufficient condition is given such that the same 
holds for holomorphic solutions at infinity. An Atkin-Swinnerton-Dyer congruence is proven for the 
coefficient series attached to Ti(7). We conclude with a consideration of asymptotics, wherein it is 
proved that many coefficient series satisfy asymptotic expressions of the form u n ~ £X n /n. Certain 
arithmetic results extend to the study of general holonomic recurrences. 

H 

^ ■ 1. Introduction and Summary of Results 

Let f(t) be holomorphic on a region Q. Then f(t) is holonomic if it satisfies a differential equation 
C with coefficients in C[i\. The set of holonomic functions is a commutative C-algebra under 
addition and pointwise (Cauchy) multiplication. Similarly, a sequence {a n } is said to holonomic 
if it satisfies a linear recurrence with polynomial coefficients. These two concepts are related: if 
fit) = ^2a n t n is holonomic, then {a n } is holonomic. The converse holds, but only in the formal 
sense. For example, {n!} is clearly holonomic. For more general theory (especially in the theory of 
S^ ■ computing with holonomic sequences and functions) see |llj . 

If the coefficients of C are polynomials in Z[i], then we say that / is holonomic over Z (and 
similarly for R, C, etc.). If f{x) is holonomic over Z[i] (or other rings), / may or may not lie in 
£T) ■ ^[[*]] (consider the functions x and e x ). 

The linear differential equations that will be of interest in this paper are the so-called Picard- 
Fuchs equations. A proper introduction can be found in [13]. Briefly, let E± be a smooth family 
of elliptic curves, parametrized over C. Fix a point to, and a cycle 7. Locally, we may identify 7 
with cycles in an open neighborhood of fibers. Now fix a holomorphic 1-form cj(t) defined over all 
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smooth fibers, and define the period of uj as 

vr(t):= fu(t). 

It can be shown that iv(t) satisfies a second order differential equation with polynomial coefficients, 
called the Picard-Fuchs equation of Et- This differential equation can be computed in full generality 
by the Griffiths-Dwork algorithm and by classical means in the case of elliptic curves. 
Consider the family of elliptic curves characterized by the Weierstrass equation 

(1) E t : y 2 + a\xy + a 3 y = x 3 + a 2 x 2 + a 4 x + a 6 
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in affine form, where aj = Oj(£) E Z[t]. We shall frequently write at for both aj(£) and aj(0), the 
reduction of at mod i. To reduce ambiguity, care will be made to differentiate between '=' and '='. 
In terms of these coefficient functions, the modular invariants g 2 and g% can be expressed as 

1202 (t) = ((a 2 + 4a 2 ) 2 - 24( ai a 3 + 2a 4 )) , 
-216gs(t) = ((a? + 4a 2 ) 3 + 36(a 2 + \a 2 )(a 1 a z + 2a 4 ) - 216(a| + 4a 6 )) . 

For convenience, we shall assume that the elliptic discriminant A(t) = g\ ~~ 27# 2 is zero at the 
origin (in any case, this can be accomplished by a simple reparametrization of t). The ^-invariant 
is written as j(t) = g 2 (t) 3 /A(t). 

It is well known that in the case of elliptic curves, the Picard-Fuchs equation is given by the 
following system of differential equations (see |14| Equation (3.3)]) 

m ± (fx\ _ 1 ( -2A'(t) 36 7 (£)\ (h 

[ ' dt \f 2 ) ~ 24A(t) V-3 52 (t)7(t) 2A'(t)J \f 2 

where 7(f) = 3gz(t)g 2 (t) — 2g 2 (t)g' 3 (t). From |13j . two independent solutions to the Picard-Fuchs 
equation at t = are 

and T(t) log(l/j(t)) +^(t), where ^(i) is holomorphic about £ = Oo Let /(£) = X^n>o n ™* n denote 
the unique holomorphic solution at t = with uq = 1. Since /(£) is a solution to a differential 
equation, the {u n } correspond to a holonomic recurrence. We have the following theorem which 
bounds the denominator of u n . 

Theorem 1.1. Let f = Xm>o Un ^ n ^ e ^ e un iQ ue holomorphic solution to the Picard-Fuchs equation 
associated to the family Et given in ([T]) and satisfying uq = 1. With d n denoting the denominator 
of u n in reduced form, d n divides the integer (8(12g2(0)) 4 ) . 

We remark that this theorem generalizes Theorem 1.5 of |13| in that Stienstra and Beukers 
specifically consider the case when a\ = 1 (mod t) and a 2 = ■ ■ ■ = oq = (mod t). A refinement 
of this theorem allows us to show that there exists a k £ Z such that f{kt) £ Z[[t]] and under 
certain additional conditions allows us to show that fit) £ Z[[£]]; that is, that the {u n } are integral, 
a result that would have been very hard to show from our recursive definition alone. To further 
strengthen this, we prove the following 'reduction' theorem, which is ultimately dependent on the 
upper bound established in Theorem 11.11 

Theorem 1.2. Fix a prime p, and let qo, ■ ■ ■ ,qe € Z[n]. Suppose that the integers {u n } satisfy the 
holonomic recurrence 

(3) (n + l) 2 u„,+i = p k °q (n)u n + p kl qi(n)u n -i + ... + p ke q e (n)u n -e, 

with uq = 1 (and we interpret U—\ = . . . = u_£ = 0) such that ki > (1 + i)ko for 1 < i < £. Then 
p |n(*D-2/(p-i))l divides u n . 

So far, our integrality results have relied upon explicit Weierstrass equations of the form (pQ), in 
which the coordinate functions a,i(t) lie in Z[t]. While this is sufficient for many applications, it 
will often be necessary to generalize this criterion to admit ai(t) £ Z(£). This is particular useful in 
conjunction with Top and Yui's work [16], in which explicit equations for higher index subgroups 
are obtained by rational maps of the parameter t. For example, the equations given for To (8) read 



Moreover, these represent a basis for solutions in any neighborhood of j(t) — 00. 
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as Et '■ y 2 = x 3 + (2 — t 2 )x 2 + x. Following the inclusion Ti(8) C T"o(8), the authors construct an 
explicit equation for Ti(8) by 

t 2 
E' t : y 2 = x 3 + (2-s 2 )x 2 +x, s=-^—. 

t z + 1 

The associated Picard-Fuchs solution, when reparametrized t \- > 4£, appears to have a power series 

with integral coefficients, where 

{u n }n=o = I 1 ; °> 0, °> 16 > 0, ~ 5l2 > °> 12864 > 0, -299008, 0, 6743040 . . .} 

yet our current extension of the Stienstra-Beukers method, Theorem II. 1\ does not permit rational 
parameters. 

Let the rational function r(t) = p(t)/q(t) G Z(i) be given in lowest terms. We define v(r) = p(0) 
and 5(r) = q(0). In general, if q(0) ^ 0, we have r(t) G Z[[i/<5(r)]]/<5(r), where the outermost 5(r) 
can be dropped if r(0) £ Z. With this new notation, we present the following generalization to 
Theorem 11.11 



Theorem 1.3. Suppose that the hypotheses of Theorem [Q] are met, with the relaxed assumption 
that the coordinates ai(t) lie in 7L(t). With f(t) = ^2 n >QU n t n and d n denoting the denominator of 
u n , it follows that d n divides the integer (8<5(A)z^(12(72) 4 ) Tl - 

We now consider Picard-Fuchs solutions at t = oo. Using ideas from the proof of Theorem 1 1.1| we 
give conditions under which the denominators of the coefficients of the Laurent series solution, in t, 
are bounded. Again under certain conditions, we can show that the coefficients of the Picard-Fuchs 
solution are integral. 

As a motivating example for our analysis of holonomic recurrences and solutions to Picard-Fuchs 
equations, we consider the solution and recurrences associated to Ti(7). This subgroup is a genus 
zero, torsion- free index 24 subgroup of PSL2(Z), and is associated with the following family of 
curves: 

(4) Et : y 2 + (1 - t - t 2 )xy + {t 2 + t 3 )y = x 3 + (t 2 + t 3 )x 2 , 
with Picard-Fuchs equation 

(5) t(t + l)(t 3 + 8t 2 + 5t- l)F" + (5t 4 + 36t 3 + 39t 2 + Bt- 1)F' + (it 3 + 21t 2 + 15* + 1)F = 0. 

Let X}n>o n n* n denote the unique holomorphic solution around t = 0. Then 

(n + l) 2 u n+ i = (2n + l) 2 u n + (13n 2 + 2)n n _! + (9n 2 - 9n + 3)n n _ 2 + (n - l) 2 n n _ 3 
(6 J 

uq = l,u\ = 1,U2 = 6, U3 = 25, U4 = 125, us = 642. 

It follows easily from Theorem 11.11 that these {u n } are integral. The family Et also admits a 
holomorphic solution about t = oo, which can be written as X^n>o v nt~ n ~ 2 with 

(n + 1) V + i = -(9n 2 + 9n + 3)v n - (Un 2 + 2)t> n _! - (2n - l) 2 v n - 2 + (n - l) 2 f n -3 
^' v = 1, vi = -3, v 2 = 12, v 3 = -59, v A = 325, v 5 = -1908. 

Again, we can show that these {v n } are integral. We can then prove that the sequence {|i> n |} 
satisfies a certain congruence relation. Using the theory of modular forms and a theorem of Verrill 
in [T7j, we obtain the following theorem. 

Theorem 1.4. Let {v n } be defined such that 

(n - l) 2 v n = (9n 2 - 27n + 21)v n _i - (13n 2 - 52n + 54)v n _ 2 + (2n - 5) 2 u n _ 3 + (n - 3) 2 v n _ 4 

with vq = 0, v\ = 1, v-2 = 3, t>3 = 12, and V4 = 59. Then for all primes p 7^ 7 and any integers m 
and r, we have 

Vmpr ~ 7pV mp r-i + (-J p 2 v mp r- 2 = (mod p r ) 
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where 7 P is the pth coefficient in the q-expansion of r](z) 5 T](7z) 3 and (§) is a Legendre symbol. 

To further study the arithmetic properties of holonomic recurrences in general, we investigate 
asymptotics of these holonomic recurrences. Consider the constant holonomic recurrence IA, given 
by 

N 

(8) U : u n+1 = ^2 kju n ^j , kj G C. 

3=0 

We define x( x ) = x N+l — k$x N — ■ ■ ■ — fcjv, typically called the characteristic polynomial of IA. Let 
Ai, . . . , X n denote the distinct roots of \- It is well known that there exist polynomials p\ t G C[n] of 
degree strictly less than the multiplicity of Aj such that u n = ^ A Xfpx^n). Thus, the asymptotics of 
IA are readily obtained, and may be simplified even further by summing over only those eigenvalues 
which have maximum absolute value. In the case of a single dominating eigenvalue A of multiplicity 
one, we obtain the particularly simple asymptotic u n ~ p\X n . 

Far less is known about the asymptotics of general holonomic recursions (i.e., when ki G C(t) in 
(jSJ)). An important subclass of such recursions are those of Poincare type, in which limn-^oo ki{n) 
is finite for all i; in that case we have the following theorems of Poincare and Perron ([61 Theorems 
8.9-8.11]): 

Theorem. Suppose that ki(n) ^ for all n G Z + . Then IA has a fundamental set of solutions 
{u 1 , . . . , u N+1 } such that limsup,^^ yfujj = |Aj|. // the eigenvalues Aj have distinct norms, then 
limn^oo u l n+1 /u l n = Aj. 

This theorem is of immediate interest due to the following proposition. 

Proposition 1.5. Let f(t) = Y2n>o u n.t n represent the power series of the solution to a Picard- 
Fuchs differential system of a family of elliptic curves. Then the holonomic recurrence of the 
coefficients of f is of Poincare type. 

Proposition 11.51 implies that the asymptotic growth of the coefficients u n in the Picard-Fuchs 
solution f(t) = Yl u nt n is approximately \u n \ ~ A™, in that for all e > we have (|A| — e) n <C 
\u n \ ^ (|A| + e) n (where '<C' denotes the asymptotic less than). In one special case, we may say 
more: 

Theorem 1.6. Let Pk(n) = a^n 2 + b^n + c^ £ Z[n] with a^ > 0, and consider the sequence {u n }, 
defined recursively by 

N 

(9) (n + l) 2 n n+ i = ^jPk{n)u n _ k . 

fc=o 

Then the characteristic polynomial \ has a unique positive real root A with multiplicity 1. Ifai~{k — 
1) + bk < for all k, then there exists an £ > such that u n < £\ n /n for n S> 0. 

Example 1.7. We illustrate the preceding with the family of elliptic curves associated to Ti(7), first 
presented in (|3|). Given the recurrence 

(n + l) 2 u n+ i = (2n + lfu n + (13n 2 + 2)u„_i + (9n 2 - 9n + 3)n„_ 2 + (n - l) 2 u„_ 3 
uq = 1, U\ = 1, U2 = 6, U3 = 25, U4 = 125, u$ = 642, 

we note that the condition a^{k — 1) + b k < is met for all k (in this case, equality holds for 
all k, which will become significant in Theorem 11.81) . Our unique positive root is A ~ 6.295897, 
a root of x 3 — 5x 2 — 8x — 1. With the notation of the proof of Theorem II. 6} let R = 21 and 
7 = A, and choose n = 10 6 . This permits j3 = 507/19, whereby we may also take £ = 0.3556365, so 
u n < 0.3556365A n /n for n»0. We will typically be interested in the infimum over all possible £ 
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in Theorem ll.6[ which can be approximated non-rigorously through an increasing sequence on n. 
This example can be strengthened with the following theorem: 

Theorem 1.8. Suppose that the conditions of Theorem \1.6\ hold, with a^{k — 1) + b^ = for all k 
and Cfc + kat > for all k. Suppose further that X is unique in absolute value, and that u n 3> for 
large n. Then u n ~ £oX n /n, for some £q G M + . 

Example 1.9. Although the hypotheses of Theorem 11.81 may seem hopelessly restrictive, we note 
that these conditions are met by the recursive sequence {u n } associated with Ti(7) (see Example 
II. 7p . Indeed, with notation as in the proof of Theorem 11.81 we have ro = 1, r\ = 15, r2 = 21, and 
r3 = 4. It follows that u n ~ £oX n /n, where A is as in Example 11.71 and 0.35561 < £q < 0.35564. 
These bounds are obtained from Example 11.71 and the equation v n > w n , coupled with the closed 
form for w n (taking M = 1000) where v n ,w n , and M are defined as in the proof of Theorem 11.81 
We note that these bounds can be taken to arbitrary precision with our current methods. 

Theorem II .81 is useful for general holonomic recurrences, however, it does not allow us to explicit 
determine the value of £q. If we restrict to studying recurrences that arise from Picard-Fuchs 
equations and hence from modular forms, one can actually deduce the asymptotic formula for u n 
with the value of £q differently. To this end, we offer the following theorem. 

Theorem 1.10. Let f{t) = ^2 n >QUnt n denote the solutions to the Picard-Fuchs equation and 
suppose the hypotheses of Theorem \1.6\ are satisfied. If the sequence {u n X~ n } is eventually positive 
and monotonically decreasing and for some constants a and b we have f{t) + alog(l — At) — )• b as 
t — > (A _1 )~ ; then u n ~ a\ n /n. 

Typically to find such constants a and b, one resorts to the theory of modular forms. As an 
application of this theorem, for the recurrence given in ([6]), we have u n ~ aX n /n where A is as in 
Examples 11.71 and 11.91 and 

7sin(47r/7) 
256vr sin 3 (6vr/7)sin 5 (27r/7) 
Analogous constants in the asymptotics associated to other congruence subgroups are given, the 
method of proving these constants is similar to that of the Ti(7) case. It seems likely that these 
constants are related to periods of certain elliptic curves. 

Similar work has been done with the index 12 subgroups (groups associated to the Beauville 
families), which offer the advantage of closed forms for the coefficients {u n }. Our analysis attempts 
to do as much as possible without the crutch of a closed form solution to the recurrence. 

Our paper is organized as follows. In Section 12. 1| we prove Theorem 11.11 and introduce the 
concept of "integral level." Sufficient conditions for integrality of certain holonomic recurrences are 
given. Theorem 11.21 is shown in Section 12.21 This is followed by Section 12.31 where we consider 
the integral levels of Ti(7) and To(12). In Section 12.41 we extend Theorem ll.il to the case when 
ai(t) £ Z(t) and prove Theorem 11 .31 We then consider the solution to the Picard-Fuchs equation at 
t = oo in Section [3l The congruence in Theorem 11.41 associated to the solution of the Picard-Fuchs 
solution at infinity is proven in Section 13.11 Finally, in Section [5] we study the asymptotics of 
holonomic recurrences, wherein Proposition 11.51 Theorem 11.61 Theorem 11.81 and Theorem 11.101 are 
proven. 

The computations for finding Picard-Fuchs differential equations (based off ([2|)) and the corre- 
sponding {u n } were done using a program the authors built in Mathematica. This program is 
available at 1 10 . 



(10) a = 3 v ' . 5 ^ w 0.3556270700876065. 



2. Near-Integrality 

Let Et be a family of elliptic curves such that A(0) = (where A is the elliptic discriminant). 
From [13] (and earlier), if Et is parametrized with coefficients in Z[i], then there exists a unique 
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holomorphic solution f(t) = Yln>o Un ^ n ^° the Picard-Fuchs equation about t = with uq = 1. 
Furthermore, f(t) is a G -function; that is, there exists a constant A such that the denominator d n 
of u n satisfies \d n \ < A n (see [5]). In what follows, we extend a result of Stienstra and Beukers [13] 
to prove an effective bound on A. 

2.1. Preliminary Bounds. We first develop some preliminary bounds on A. To do this, we will 
need the following lemma. 



Lemma 2.1. For integral k > 0, we have 

12 



3kl -5/12V-V12 



fc / V k 

Proof. Up to units, we expand the product as 12 k /(kl) 2 FL=o(^ + 12j)(l + 12j) . For prime p, let 
n p and dp denote the exponent of p in the numerator and denominator, respectively, in the above. 
Then n 2 = 2k, while Legendre's formula implies that 

oo 

d 2 = 2([k/2\ + [k/4\ + [k/8\ +...)< J2 k / 2i = 2k ■ 

i=0 

Similarly, one may show that n% = k > d%. Now let p > 3 be prime. We have p | (5 4- 12j) if and 
only if j = — 5(12~ 1 ) (mod p), and a similar formula holds for (1 + 12j). Running through j up 
to k — 1, we have n p > 2(|_(fc — l)/pj + |_(& ~~ 1)/.P 2 J +•••)• Furthermore, if k = (mod p), then 
n p {k) = n p {k + 1) as functions of fc, as neither (5 + 12/c) nor (1 + 12/c) will be divisible by p. It 
follows that n p > 2([k/p\ + [k/p 2 \ +•••). As the right hand sum is d p , we conclude that n p > d p 
for all primes p, whence our product is integral. □ 

Proof of Theorem \l.l\ Throughout, let F(a,b,c;z) denote the hypergeometric function 2 Fi. From 
[13] . the function 

gives a solution to the Picard-Fuchs differential equation for E% in a neighborhood of j(t) = oo, 
which holds at t = since 52(0) 7^ and A(0) = 0. By the preceding remarks, 12^2 (t) G Z[£], so 
we have 12^2 (t)/( 12^2(0)) G Z[i/(12g2(0))], using the fact that our constant term is 1. We recall 
that (1 + 82)- 1 / 4 G Z[[z\], hence (12 52 (*)/125 2 (0))" 1/4 G Z[t/ (96^(0))] under composition!! (Note 
that all convergence properties are satisfied). 

Next, note that l/j(t) = 12 3 A(t)/(12g 2 (t)) 3 lies in l2 3 Z[[t/(l2g 2 (0)) 4 ]}. To see this, we observe 
that (12g 2 (t))- 1 G Z[[t/(12g 2 (0))]]/(12g 2 (0)), hence (12g 2 (t))- 3 G Z[[i/(12 5 2(0))]]/(12<7 2 (0)) 3 . One 
may show that A(i) G Z[t], whereby A(t) (12 g 2 (t))- 3 G Z[[t/(12 52 (0))]]/(12 ff2 (0)) 3 as well. Since 
A(0) = 0, our constant term is and we have l/j(t) = 12 3 A(t)/{12g 2 (t)) 3 G 12 3 Z[[i/(12 52 (0)) 4 ]]. 
By Lemma EH it then follows by composition that F(^, i, 1; l/j(t)) G Z[[t/(125 2 (0)) 4 ]]. Multi- 
plication then gives f(t) G Z[[i/(8 • (1252(G)) 4 )]], as claimed. This completes the proof of Theorem 

rrn '" n 



Remark 2.2. As 12g 2 (t) lies in Z[t], Theorem 11.11 will never suffice to prove integrality; we may 
prove, at most, that f(8t) G Z[[t]], under the assumption 12^(0) = dbl. In [13] on the other hand, 
Stienstra and Beukers treat a less general case, in which a\ = 1 and a 2 = . . . = oq = 0. With these 
assumptions integrality can be shown directly, upon careful examination of the expansion 

(12) 12g 2 (t) = (a^ 1 ■ (1 + Aa 2 a- 2 y l l 2 ■ (1 - 24( ai a 3 + 2a 4 )(a 2 + Aa^" 2 )- 1 ^ 



The function (1 — 8z) 1 ' 4 is the generating function for the sequence {a n } given by a n := (2"/n!) FIfc=o (^ "f" -0- 
Integrality of this sequence follows in the manner of Lemma |2. II 
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Herein, we have simply traded strength for generality. In any case, we have the following: 

Corollary 2.3. Let f = Yln>o Un ^ n as before. Then a finite number of primes divide elements of 
the denominator set {d n }. Moreover, there exists a k £ Z such that f(kt) G Z[[i]]. 

Corollary 2.4. With the notation of Theorem \1.1\ suppose that g 2 (t)/g 2 (0) is a fourth power in 
Z[[t]] and that 12g 2 (0) = ±1. Then the Picard-Fuchs solution f(t) lies in Z[[i]]. 

Proof Our additional assumption implies that a(t) = {g 2 {t)/g 2 (0))~ 1 ' 4 € Z[[i]], using the fact that 
a(0) = 1. For the hypergeometric term in equation (fTTj) . note that l/j(t) = 12 3 A(t)/(12<72(i)) 3 lies 
in 12 3 Z[[i]], whereby Lemma O implies F(^, ±, 1; l/j(t)) £ Z[[i\] as well. □ 

Remark 2.5. This last corollary is the essence of Stienstra and Beukers integrality result in Theorem 
1.5 of |13j . wherein the decomposition of (|12p is used to show that our 'fourth power' hypothesis is 
met. 

Example 2.6. Consider the family of elliptic curves parametrized by 

E t : y 2 + 5xy + (t 2 + l)y = x 3 + (t - 3)x 2 + x + t 3 + 1 . 

We compute 12g 2 (t) = 1 + 104t - 104t 2 and A(0) = 0. With a(t) = I3t - 13, we may write 
(12g 2 {i))~ l ' A = (1 — 8ta(i)) -1 ' 4 . As in the footnote of the proof of Theorem 1 1.1^ this implies that 
12<72 (t) has a fourth root in Z[[t]]. By Corollary 12.41 it then follows that the Picard-Fuchs solution 
f(t) = Y2 n >o u nt n nes i n ^[Wl- The (integral) coefficient series {u n } satisfies a 12-term holonomic 
recurrence and begins with 

1, -86, 26856, -10713740, 4757138560, -2243385809196, 1099636896720096... 

This example is taken only to illustrate that the requirements a\ = 1 and 02 = • • • = «6 = from 
[13j are not strictly needed in elementary proofs of integrality. In Example 12.81 a more general 
example is introduced. 

A parametrization of the family Et is said to be reduced if a,i{t) € Z[ni] for all % implies that 
n = ±1. Clearly, every family is equivalent to a reduced family with A(0) = 0, under linear 
transformation. By Corollary 12.31 every reduced family admits a reparametrization 1 1— > kt, k 6 Z, 
such that f(kt) G Z[[t]]. 

Definition 2.7. If E has a reduced parametrization and k is minimal (in absolute value) such that 
f(kt) £ Z[[t]], then we say that \k\ is the integral level of E. 

One may extend this with the additional notion of fractional integral level, corresponding to the 
cases in which u n /k n £ Z for all n. In any case, these will be well-defined by our previous work, 
and should be viewed relative to the reduced representative of an elliptic curve equation. Similar 
adjectives will be used to describe the associated recurrences in the coefficients u n of f. 

Example 2.8. Consider the (reduced) family of curves parametrized by 

E t : y 2 + 2xy = x 3 + x 2 + tx + t 

We calculate A(i) = -16t(4t 2 - 13* + 32) and g 2 {t) = 16/3-4*. By Theorem[LTJ we have d n \ 2 27n , 
or equivalently, /(2 27 t) G Z[[i]]. The coefficient recurrence associated to f(2 27 t) is given by 

(n + l) 2 n n+ i = 2 20 (36n 2 + 68n + 9)u n - 2 44 (76n 2 + 104n - 123)w n _i - 2 71 (16n 2 - 48n + 35)n n _ 2 , 

u = 1, m = 2 20 3 2 , u 2 = 2 38 105 

hence {u n } is an integer sequence, despite division by (n + l) 2 at every iteration. To contrast this 
with the 'average' behavior of similar holonomic sequences, see Zagier's work |20| . We stress that 
this example does not imply that the integral level of E is 2; however, given our previous work, it 
does show that the level of E divides 2 . On the other hand, one may check (by inspection) that 
the level of E is at least 2 9 . 
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Remark 2.9. The Weierstrass form of an elliptic curve is not unique up to isomorphism, but one 
may show that the j-invariant is, and that the modular invariant g<i is unique up to multiplication 
by a fourth power |12j . Expanding this to families of elliptic curves, the formula in equation 
(jlip implies that the solution to the Picard-Fuchs equation is unique up to linear transformation 
and multiplication by a polynomial in Z[i]. To avoid this technical matter in what follows, explicit 
equations will be given for our parametrized families of curves. Note, in particular, that any scaling 
of the argument of / (e.g. during the integral level analysis), corresponds to a trivial action on the 
isomorphism classes of families of elliptic curves. 

2.2. Further Bounds. In this section we give further bounds on the integral level of Ef. We first 
give the proof of Theorem 11.21 

Proof of Theorem \1.2l For integer n, let p(n) denote the maximal power of p dividing n. By ([3]), 
we have 

2p(n + 1) + p{u n+ i) > min (p(u n ) + fa,p(u n -i) + fa, . . . , p(u n _ e ) + fa) 

(13) > fa + min (p(u n ), p{u n -i) + fa-fa,..., p{u n ^i) + fa - fa) 

As a special case, suppose that the minimum in (I13D is always obtained by the first term, p(u n ). 
Then 

p(u n +i) >fa~ 2p(n + 1) + p(u n ) > 2fa - 2p(n + 1) - p(n) + /o(« n _i) > . . . 

n 
> fa(n + 1) - 2^2 p{n + 1 - i) + p(u ) > k (n + 1) - 2p({n + 1)\) 
i=0 

wherein we have used the 'logarithmic' property of p and the fact that p(uq) = 0. By Legendre's 
formula, we have p((n+ 1)!) < (n + l)/(p — 1), hence p{u n+ \) > (n + l)(&o — 2/(P — I))- Integrality 
of p then implies that p{u n+ i) > \(n + l)(fa — 2/(p — 1))]. This concludes the special case. Next, 
suppose that the minimum in the first iteration of (1131) is the term p(u n -i) + fa — fa, instead of 
p(u n ). Then 

p{u n +i) > fa- 2p(n + 1) + (p(u n -i) + fa- fa) > (i + l)fa - 2p(n + 1) + p{u n -i) 

>{i + l)fa - 2{p(n + l) + ... + p(n + l-i))+ p(u n -i) 

Note that this is precisely the inequality we would have derived after (i + 1) iterations of the above 
induction. In general, this deviation from the special case has the same effect (none) regardless of 
when it is implemented in the algorithm. As the index of Uj falls below £, note that the recursive 
definition in ([3|) will begin to omit summands, beginning with j r e qg(n)u n -g, hence the algorithm 
will be forced to conclude along the lines of the special case (instead of jumping over no to negative 
indices) . This completes the proof of Theorem 11.21 □ 

While the previous section gives a theoretical upper bound for the integral level of a family of 
elliptic curves, Example 12.81 suggests that our current bounds are far from optimal. 

Corollary 2.10. With {u n } as in ©, let s = \(fa - 2/(p - 1))] . Then v n := p~ ns u n G Z, and the 
sequence {v n } satisfies the recurrence relation 

(14) {n + l)V+i = p k °- s q (n)v n + p k ^ 2s q x {n)v n ^ + ... + p ke " {e+1)s qe(n)v n . e . 

In particular, f{p~ s t) G Z[[t]]. 

Example 2.11. We continue ExampleES Taking s = 18, we find that f(2 27 ~ 18 t) = f(2 9 t) G Z[[t]], 
with coefficients {u n } satisfying the integral holonomic recurrence 

(n + l) 2 u n+1 = 2 2 (36n 2 + 68n + 9)u n - 2 8 (76n 2 + 104n - 123)n n _i - 2 17 (16n 2 - 48n + 35)-u n _ 2 , 
uq = 1, u\ = 36, U2 = 420. 
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It follows that E has integral level exactly 2 9 , where the lower bound stems from Example 12.81 
While our algorithm here gives definite results, note that this does reflect the general case. 

2.3. Some Applications. In this section, we apply our analysis to several curves associated to cer- 
tain genus zero congruence subgroups of SL2 (Z) . By Remarks 12.21 and 12.51 the following proposition 
is immediate: 

Proposition 2.12 (cf. (1.5) in [13]). Let Et be a family of elliptic curves, parametrized as in (fT]). 
Assume A(0) = and let f(t) = Yln>o Un ^ n ^ e ^ e un iQ ue holomorphic solution to the Picard-Fuchs 
equation with no = 1. Suppose ([T]) reduces modulo t to y 2 ± xy = x 3 . Then u n G Z for all n > 0. 

We now apply the above proposition to prove the integrality of two sequences related to the 
congruence subgroups Ti(7) and To(12), respectively. In the following, we shall always take f(t) = 
Y^ n >o u nt n , with u = 1. 

Example 2.13. Consider the family of elliptic curves associated to the congruence subgroup Ti(7) 
given in Q with the associated recurrence in (|rj]). Observe that (JU) satisfies all the conditions of 
Proposition 12.121 Then m„£Z for all n > and Et has integral level 1. 



Example 2.14. Next, we consider the family of elliptic curves associated to To(12). Again from |16| . 
the parametrized Weierstrass equations for this family of elliptic curves is given by 

(15) Et : y 2 + (t 2 + l)xy - t 2 (t 2 - l)y = x 3 - t 2 (t 2 - l)x 2 . 

Accordingly, the Picard-Fuchs equation is 

t(9t 4 - 10t 2 + l)F" + (45£ 4 - 30£ 2 + \)F' + 12£(3£ 2 - l)F = 0. 

The recurrence for this example is 

(n + l) 2 n n+ i = (10n 2 + 2)« n _i - 9(n - l) 2 n n _ 3 
(16) 

no = 1, u\ = 0, U2 = 3, n 3 = 0, m = 15. 

Once again, the hypotheses of Proposition 12.121 are met, hence u n £ Z for all n > 0. 

Remark 2.15. Another way to prove u n G Z is to observe that U2 n +i = for all n and ui n = 
Sfc=o (fc) ( k ) ' unfortunately, closed forms of this type (sums of products of binomial coefficients) 
have not been found for other recurrences derived from congruence subgroups of level greater than 
120 

2.4. Integral Analysis with Rational Functions. In the previous sections, we assumed the 
Oj(i) € Z[i]. Now we generalize the work in Section [2.11 to when di(t) G Z(i). 

Proof of Theorem \1.3X We proceed along the lines of the proof of Theorem 1 1.11 From the notational 
remarks preceding the statement of Theorem II. 3| it follows that 

12<fc(t)/(120j(O)) G Z[[t/6(12g 2 /(12g 2 (0)))]] C Z[[t/(12g 2 (0)5(12g 2 ))}} = Z[[t/u(12g 2 )]}. 

As before, we obtain (12 5 2(t)/(12g 2 (0))) _1/4 G Z[[t/(8v(12g 2 ))]]. 
Next, note that (12g 2 (t))- 3 G Z[[t/^(12o 2 )]]/^(12^ 2 (0)) 3 , hence 

(17) ^ - 12 A(t)(12g 2 (t)) G 12 -^-^ . 

As A(0) = 0, the external denominator in (117j) can be absorbed into the brackets, and we 
have (by Lemma 12. ip that the hypergeometric term of the Picard-Fuchs solution f(t) lies in 
Z[[t/(5{A)v(12g 2 )*)]\. Upon multiplication we derive f{t) G Z[[t/fc]], where k = 85(A)u(12g 2 ) 4 . 
This completes the proof of of Theorem 11.31 □ 



This technique appears to work in the case of To (12) for reasons described in [16] - that To(12) may be obtained 
as a double cover of the congruence subgroup To (6) (which is of level 12). 
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The bound given in the statement of Theorem 11.31 may be at times quite large. In any case, 
Theorem 11.31 does provide us with an existence result, along the lines of Corollary 12.31 Moreover, 
we remind the reader that the results of Section 12.21 apply equally well in our new situation, which 
will immediately reduce our bound on the integral level to a manageable, if not optimal, magnitude. 

Additionally, we note that Theorem 11.31 implies Theorem II .11 Indeed, when the coordinates ai(t) 
are relegated to Z[i], we have v(\2g2) = 12^2(0) and 5(A) = 1, and we may recover our earlier 
bound. Our work in Theorem 11.31 also applies to Corollary 12.41 which we now formalize as 

Corollary 2.16. With the notation of Theorem \1.3l suppose that ff2(£)/<72(0) G Z(i) is a fourth 
power in Z[[i]]. Then the Picard-Fuchs solution f(t) lies in Z[[i/(<5(A)z/(12<72) 4 )]]- 

Along the lines of Remark 12.51 we now offer a sharper version of the Stienstra-Beukers bound 
given in [13]. While technically weaker than Corollary 12. 161 its hypotheses are verifiable by inspec- 
tion alone. 

Corollary 2.17. Suppose that the family of elliptic curves Et meets the hypotheses of Theorem 
\1.S\ with the additional assumptions that 

(i) E t reduces mod t to y 2 + xy = x 3 
(ii) 5 (at) = ±1 for % = 1, . . . , 6. 

Then the integral level of E is 1. 

Proof. One may compute that 5(A) divides 5(ai) 6 5(a2) 3 5(a3) 4 5(a4 : ) 3 5(aQ) 2 in general, hence (ii) 
implies that 5(A) = 1. By condition (i), we find 1/(12*72) divides ai(0) 4 5(a2) 2 5(a3)5(a4 : ), which 
again is 1 by hypothesis. Moreover, (i) implies that Corollary 12.161 applies, which now gives our 
desired result. □ 

Example 2.18. In the introduction, we mentioned a family of elliptic curves E' t , associated to the 
congruence subgroup Ti(8). Explicitly, we have 



E' t : y 2 = x 3 + (2 - s 2 )x 2 + x 



t 2 



t 2 + l 



We compute 5(A) = 1 and 1/(1252) = 16, hence Theorem 11.31 implies that f(2 t) G Z[[i]]. The 
coefficients u n of f(2 19 t) satisfy 

(n + l) 2 u n +i = -2 39 (2n 2 - 3n + l)u n _i - 2 74 (23n 2 - 118n + 143)u n _ 3 

- 2 113 (7n 2 - 61n + 130)u n _ 5 - 2 150 • 3(n 2 - 12n + 35)u n __ 7 , 

and so an appeal to Corollary 12. 101 with ko = 18 (whence s = 16) implies that f(8t) G Z[[i]]. On the 
other hand, us = 201/1024, which implies that a scaling t 1— >• 4t is necessary to ensure integrality. 
Thus, the integral level of E' t is either 4 or 8. 

A bit more can be said for this example. Based on the recurrence relation on {u n }, we deduce 
that / is an even function. Setting v n = U2 n , we derive the following integral recurrence for 
8 2n u 2n = ^ n v n : 



(n + 1) V+i = 2 b (4n 2 + n)v n -2 w \2?>n 2 - 36n + 12)v n _ 1 

-3 • 



+ 2 16 (14n 2 - 47n + 38)f n _ 2 - 2 22 (3n 2 - 15n + 18)u n _ 



This sequence admits a reduction (take kg = 5 in Corollary I2.10p . whence we deduce that 8 n v n G Z. 
This implies that 4 n n n G Z for all n, whereby /(4i) G Z[[t]] and we have shown that the integral 
level of E' t is precisely 4. Indeed, we have the stronger result that f(2\/2t) G Z[[t]], which - while 
not an integer scaling - remains significant for number theoretic reasons. 
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3. Picard-Fuchs Solutions at Infinity 

Let Et denote a family of parametrized elliptic curves, with j-invariant j(t) = g 2 (t) 3 /A(t). Thus 
far we have considered only the holomorphic solution to the Picard-Fuchs equations at t = 0. 
This naturally extends to power series solutions centered at i G Z (where A(i) = 0) by affine 
transformation, which will preserve integrality in our power series coefficients (as a composition of 
functions in Z [fill). To extend these results to the case i = oo, we recall that the function 



(18) f{t) 



(12g 2 (t))- 1 / 4 F(—,-,l;4r) 
v yi\ )) ^12' 12' j(t)J 



of Theorem 11.11 gives a solution to the Picard-Fuchs equation for E% in any neighborhood of j = oo 
(although, in general, this will not give uq = 1 as we have scaled differentially here compared to 
(fTT]) ). This, of course, does not require A(i) = 0; instead, we take i = oo and make throughout 
this section the assumption that j(oo) = oo. 

Theorem 3.1. Consider the parametrized family of elliptic curves given by 

E t : y 2 + ai(t)xy + a 3 (t)y = x 3 + a 2 (t)x 2 + a 4 (t)x + a G (t) , 

where a% G Z[i]. Let a^ denote the leading coefficient ofl2g2(t), and suppose that j'(oo) = oo. If 
4 | degg 2 , then the Picard-Fuchs differential equation for Et admits a holomorphic solution of the 
form 

f(t) = Y,U n t- n - {Acg92)/i 
n>0 

in a neighborhood of t = oo, in which we may take Uq = 1. Let d n denote the denominator of u n in 
reduced form. Then d n divides (8ajj) n . If V2g 2 (l/t)t deg92 /a^ is a fourth power in Z[[i]], then we 
have the stronger result d n \ aj^. 

Proof. Let s = 1/i. Then [I2g 2 {\/ s)s deg92 )~ l G Z[[s/a N ]]/a N , and hence (12g 2 (l/s)/a N )- 1 G 
s deg92 Z[[s/aj\r]] by the remarks preceding the statement of Theorem 11.31 Since 4 | degg 2 , it 
follows that {12g 2 (l/s)/a N )- 1 ^ G s^ dcg92 ^ 4 Z[[s/(8a N )}} as in Theorems O and O Moreover, 
if (12g 2 (l/s)s deg92 /ajv) is a fourth power in Z[[s]], then the stronger result (12g 2 (l / s) / aj^)^ 1 ' 4 G 
s^ deg92 ^ 4 1[[s/a N }} also holds. 

Next, we consider l/j(l/s) = !2 3 A(l/s)(l2g 2 (l/s))~ 3 . Within this factorization, we have 
(12g 2 (l/s))- 3 G s 3deg92 Z[[s/a N ]]/a N and A(l/s) G s" degA Z[[s]]. As j(oo) = oo, we have that 
l/j(l/s) ->■ as s — > 0. Thus, in the product A(l/s)(12g 2 (l/s)~ 3 ), we have 3degg 2 - deg A > 0. 
Using this, we conclude that 

— l — G l2 3 s 3dcg92 ~ dcgA Z[[s/a%}}. 

In other words, the external denominators arising from (12g 2 (l/s))~ 1 may be brought within the 
brackets. Upon multiplication we derive F(^, ^, 1; -/h ■ , ) G Z[[s/a^]], hence the Picard-Fuchs 
solution /(i) given by 



\ a N J ' V 12 ' 12 ' 'i(0, 

lies in s^ deg92 ^ 4 Z[[s/(8a 4 N )]}. Furthermore, /(*) G s^ deg92 ^ 4 Z[[s/a%\] when a N /(l2g 2 (t)) is a fourth 
power in Z[[s]]. In either case, we may write f(t) = ^2 n>0 u n t~ n ~^ dog92 '' 41 , in which we now 
claim that uq = 1. Our assumption j(oo) = oo implies that lim s _ > o^ ? (r2' 12 ' ^> 7(T7sT ) = ^> nence 

uq = liuis-^os^ S92 '' (12g 2 (l/s)/aN)~ 1 ' 4 - As such, we may take no = 1 in the expansion of 
/(i) and need not worry about the effects of rescaling on integrality. This completes the proof of 
Theorem 13.11 □ 
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Several of our previous results have relied on the hypotheses that (<72(*)/S2(0)) or 12^2 (i) / '&N 
admit fourth roots in Z[[£]] (see Corollaries 12.41 and 12. 161 Theorem 13. ip . Unfortunately, the authors 
are unaware of any property or condition that is necessary and sufficient to show this result. As 
such we have only addressed two special cases thus far: 

(i) that which follows from ai = 1 and 02 = ■ ■ ■ = a& = 0, or 
(ii) that which follows from (1 + St)- 1 ^ G Z[[t]]. 

The following proposition, taken from a recent paper of Heninger, Rains, and Sloane [8], can be 
used to extend these techniques: 

Proposition 3.2 (cf. Theorem 1 in |8j). Fix an integer n, and define \x n = wJT_i n p, where the 

product ranges over the primes dividing n. Let c(t) = ^Ckt € Z[[£]], with cq = 1. Then c and c 
(mod Hn) admit nth roots in Z[[£]] simultaneously. 

From this we derive the following trivial corollary. 

Corollary 3.3. Let pit) £ Z[i] with p(0) = 1, and suppose that the coefficients of p reduce modulo 
8 to the fourth power of a polynomial in Z[i]. Then p(t) admits a fourth root in Z[[i]]. 

Example 3.4. We recall from Example 12.131 that the family of elliptic curves 

E t : y 2 + (1 - t - t 2 )xy + (t 2 + t 3 )y = x 3 + (t 2 + t 3 )x 2 

corresponding to Ii(7) has integral level 1. Based on the calculation 

(t 8 + 12t 7 + 42£ 6 + 56t 5 + 35£ 4 - 14t 2 - At + l) 3 



m 



UH 7 (t + l) 7 (t 3 + 8t 2 + ht - 1) 



and the fact that deg<?2 = 8, it follows that Et admits a holomorphic expansion about t = 00, 
which fits the form f(t) = ^2 n >QUnt~ n ~ 2 by Theorem 13.11 In the notation of Theorem 13.11 we 
have on = 1, hence the integral level of Et at infinity divides 8. To strengthen this, we appeal to 
Corollary 13.31 we have 

12t 8 g 2 (l/t)/a N = t 8 - At 7 - Ut 6 + 35t 4 + 56t 3 + 42t 2 + 12* + 1 = (t 2 + 1 + l) 4 (mod 8), 

whence Corollary 13.31 implies that 12t 8 g2(l/t)/aN is a fourth power in Z[[i]]. By this, Theorem 13. II 
now shows that f(t) E Z[[i]], so Et has integral level 1 at infinity. The associated recurrence of the 
holonomic recurrence {u n } is given by 

(n + l) 2 u n+ i = -(9n 2 + 9n + 3)u n - (13n 2 + 2)u n _i - (2n - l) 2 u n - 2 + (n - l) 2 u n _ 3 
(19) 

no = 1, U\ = —3, U2 = 12, U3 = —59. 

This sequence (up to sign) will be studied in the next section. One should note the large similarity 
between the {u n } recursion and the recursion given in equation ([6]). This is no accident; whenever 
the Picard-Fuchs equation for a family Et admits singular expansions fo(t) = ^ u n t n and foo(t) = 
^2v n t~ n ~( dQg92 '' 4 about t = and t = 00, respectively, the holonomic recursions {u n } and {v n } 
are related by the following proposition: 

Proposition 3.5. Let{u n } and {v n } be as above. Lf{u n } satisfies the finite holonomic recursive re- 
lation Y^k w k(n)u n -k = for some Wk G Z[n], then {v n } satisfies the recursion ^ fc Wk(— n)v n+ k+s = 
0, where S > is the minimal integer such that v$ ^ 0. Under these recursions, the sequences {u n } 
and {v n } are uniquely defined by the conditions Uq = v§ = 1, and u n = v m = for n < and 
m < 5. 

Proof. Let p{t)f{t) + q(t)f'(t) + r(t)f"(t) = denote the associated Picard-Fuchs equation, where 
p(t) = J2pit l , q(t) = ^qit 1 , and r(t) = ^ r^f 1 are polynomials in Z[i]. From the solution /o we 
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derive 

^2 ( pk ~ kqk + 1 + r k+2(k + l)k + n(q k+1 - r k+2 (2k + 1)) + n 2 r k+2 ) u n _ k = , 
k 
seen by equating the t n coefficient to and regrouping. Similarly, the solution f^ implies 

^2 ( pk ~ kqk + 1 + r k+2{k + l)k + n((2k + l)r k+2 ~ Qk+i) + n 2 r k+2 ) v n+& +k = , 

k 

upon consideration of the t~ n ~ term. (For simplicity of notation, one should take these sums as 
infinite with finite support.) Our first result follows by taking w k {n) = p k — kq k+ i + r k+ 2(k + l)k + 
n((2k+l)r k+ 2—Qk+i)+n 2 r k+ 2- The initial conditions are trivial; merely stated for completeness. □ 

3.1. A Congruence Associated to Ti(7). We prove Theorem ll. 41 which gives an Atkin-Swinnerton- 
Dyer congruence for the sequence (up to sign) given in (|19p . In this section, we let q = e 2mT where 
Im(r) > 0. Note that 

oo 
/( T ) = q - q Z + 2q 4 + 2q 5 - 3q & + q 7 + 3q 8 + • • • = q J](l " <fT" 

n=l 

where c n = 2,0,1,-2 if n = 0, ±1,±2, ±3 (mod 7), respectively, is a weight 1 modular form for 
Ti(7). The function 

oo 
t(r) = q - 3q 2 + 5q 3 - 6q 4 + 7q 5 - 7q 6 + 3g 7 + 4q 8 + ■ ■ ■ = q J^[(l - q n ) dn 

n=l 

where d n = 0, 3, —2, — 1 if n = 0, ±1, ±2, ±3 (mod 7), respectively, is a Hauptmodul for 1^(7) (cf. 
(4.23) in [7jjj and hence is weight modular function. In the following we shall sometimes abuse 
notation and write f(q) and t(q) instead of /(r) and £(r). 

Define {v n } such that /(t) = J2 n >o Vn ^( T ) n - Computation yields that the inverse series q(t) of 

t(q) is q(t) = t + 3t 2 + 13t 3 + 66t 4 + 365£ 5 + 2128t 6 + 12859t 7 + 79745£ 8 H and hence from how 

f(q) is defined above, we have 

f(t) =t + 3t 2 + 12t 3 + 59t 4 + 325t 5 + 1908t 6 + 11655t 7 + 73155t 8 + • • • . 

The following proposition shows that these coefficients agree (up to sign) with the sequence from 

csD. 



Proposition 3.6. The {v n } are such that 

(n - lfv n = (9n 2 - 27n + 21)u„_i - (13n 2 - 52n + 54)w n _ 2 + (2n - 5) V„ 3 + (n - 3) 2 v n - 4 

with vo = 0, v\ = 1, V2 = 3, V3 = 12, and v& = 59. 

Proof. Define the two differential operators D t = t-^ and D q = q-j-. Let G\ = (D q t)/t, G2 = 
(D q f)/f, 

D q G\ - 2GiG 2 , , D q G 2 — G 2 
PiW = q2 ' and Paw = q2 • 

Then by Theorem 1 of [18 i, D 2 f + Pl D t f + p 2 f = t 2 f" + tfa. + 1)/' + P%f = 0. Since t is a 
Hauptmodul, p\ and P2 are both rational functions of t. Observe that 

2t 4 + 4t 3 -9t + 2 , ,. (t 2 - t + l){t 2 + 2t - 1) 

Pl(*) = u T\77A I CJ 9 ol I in and P2W ~ 



(t-l)(i 3 + 5i 2 -8t + l) ^ zw (t-l)(t 3 + 5t 2 -8t + l)' 



The Hauptmodul that Elkies takes is (in his notation) y 2 z/x 3 where x, y, and z are defined as in (4.4) of his 
paper. Since the Hauptmodul is the generator of the function field of Xi(7), we have taken our Hauptmodul to be 
-x 3 /(y 2 z). 
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It follows that 

t 2 (t - 1)(£ 3 + 5t 2 - 8t + l)f" + t(3t 4 + 8t 3 - 13£ 2 + 1)/' + (t 2 -t + l)(t 2 + 2t - 1)/ = 

and hence we have the desired recurrence. This completes the proof of Proposition 13.61 □ 

Remark 3.7. One can also prove the above proposition by methods similar to the discussion on 
Pages 58-60 of [2] and the values of £(r) at the inequivalent cusps of Ti(7) from Table 6]. 

To prove Theorem 11.41 we use the following theorem of Verrill: 

Theorem 3.8 (cf. Theorem 1.1 in [IT]). Let T be a level N congruence subgroup for SLi2(Z). Let 
t(r) be a weight modular function for V. Let /(r) be a weight k modular form for T, and let g(r) 
be a weight k + 2 modular form for T which has an Euler product expansion. Suppose that we can 
write 

fir) = £ b n t(r) n and g(r) = £ ln q n 

n>0 n>0 

where q = e 2mT . Suppose that for some integers M and a^, d \ M, we have 

^ d\M 

then for a prime p \ NM and integers m, r, we have 

b mp r - J P b mp r-i + e p p k+l b mp r-2 = (mod p r ) 

where e v is the character of g. In particular ifb\ = 1, then b p = j p (mod p). 

Proof of Theorem \1.4\ In the notation of Theorem 13.81 take g(r) = ry(r) 3 ry(7r) 3 where tj{t) is the 
Dedekind eta function. Since the weight 3 cusps forms for Ti(7) (of which g(r) is a member) 
has dimension 1, any modular form in it is a Hecke eigenform and hence has an Euler product 
expansion. Also note that in this case, the character of g(r) is the real nontrivial character mod 7 
and hence e p = (£). A finite computation yields that 

,/ x Qdt 
f{T) -ld- q =9{T) 



1pV mp r-i + (■=) V 2 v mp r-2 = (mod p r \ 



and hence we have 

for all integers m and r and p / 7 with 7 P the pth coefficient in the q-expansion of t/(t) 3 ?7(7t) 3 . 
This completes the proof of Theorem 11.41 □ 

4. Asymptotic Behavior 

We now investigate asymptotics of {u n } associated to Picard-Fuchs equations and holonomic 
recurrences in general. To do this, we first give the proof of Proposition 11.51 

Proof of Proposition \ 1.51 Let p = ^2pit l , Q = Y1 1i^-> an d r = ]C r $ % ^ ^-W denote the coefficients 
in the Picard-Fuchs differential equation p(t)f(t) + q(t)f'(t) + r(t)f"(t) = 0. Equating coefficients 
for t n gives 

deg(p) deg(g) deg(r) 

(20) ^2 PiU n -i + ^2 Qi(n + l- i)u n+1 ^i + ^ n(n + 1 - i)(n + 2 - i)u n+2 -i = 

i=0 i=0 i=0 

One may show that each zero of A(i) is a zero of r(t), hence ro = and the u n +2 term in (|20p 
is thus not present in the holonomic recursion for {u n }. We note that no u n -i coefficient in (|20p 
has degree (in n) exceeding 2. It now suffices to show r\ ^ 0; it would then follow that the degree 
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of the coefficient of u n +i (the leading term) will not be strictly exceeded, hence {u n } will be of 
Poincare type. 

Suppose t\ = 0, so r has a root of order 2 at t = 0. As the Picard-Fuchs equation is an ordinary 
differential equation of Fuchsian type, we must have go = (else the differential equation would 
admit a non-regular singular point at t = 0). Taking p, q, and r coprime without loss of generality, 
we then have po ^ 0. Equating the constant term in (j20|) gives po^o + QoUi + vqU2 = 0. Thus 
PoUo = 0, and it follows that uq = 0. This contradicts the assumption that uq = 1 in the canonical 
holomorphic solution to the Picard-Fuchs equations. Thus r\ ^ 0, and we have shown that {u n } is 
of Poincare type. This completes the proof of Proposition 11.51 □ 

Proof of Theorem \1.6l We have x( x ) = x +1 — aox N — ... — a^^\x — a at. As not all aj are 0, we 
may factor x as x( x ) = x m p(x), in which p(Q) is strictly negative. Thus, for some e > 0, we have 
p{e) < 0, hence x{ £ ) < 0. As % is asymptotically positive, let A denote the minimal positive root 
of x- For r = 0, . . . , N we have 

(N + l)\ ( x N + i_y ak (N-k)l (N + l 




X r (N + l-r)\ \ ^ (N + l)l (N-k 

As all derivatives of \ are non-negative on [A, oo), it follows that x ls strictly increasing on that 
interval and that A is the unique positive root of x (with multiplicity 1). Next, define v n := nu n /X n , 
and let Pk{n) = a^ ■ (n + l)(n — k) + r^(n), so that deg(r^) < 2. Then (|9|) implies that the sequence 
{v n } satisfies 



N , N N N 



V n +X 



p k (n)v n - k sr^ a kV n -k , v^ r k(n)v r , 



E Pk\n)V n -k _ ST^ UkVn-k V^ 

\ k + 1 {n-k){n + l)~ ^ X k+1 ^ \ k {n-k)(n + l) 



We note that the hypothesis at{k — 1) + bk < for all k implies that the polynomials rk{n) have a 
non-positive linear term. Thus, take R > such that R > r^(re) for all n > 0. We have 

^ a fc ^ Rmax.(v n ,... ,v n - N ) R{N + l)m n 

^<LAW max K-.^)+L A*+i(n-fc)(n + l) - m " + 7 (n-A0(n + l)' 

where m^ = max(f n , . . . ,v h _n) and 7 = min(A, . . . , A JV+1 ). Note that the reduction of the first 
sum to m n follows from the (algebraic) definition of A. Take /3 > such that t> n +i < m n + (3m n /n 2 , 
and choose £ > such that v n , . . . , v n —N satisfy «j < £e~ 2 ^' n for some n > max(3, /3/2) (allowed to 
be arbitrarily large). We claim that t> n +i < £e~ 2 ^'' n+1 '. To see this, it suffices to show £e~ 2 ^' n + 
£pe~ 2 P/ n /n 2 < ^ e -2/3/("+i) 5 or equivalents, 

(21) e -2/3/(n+l) _ e -2P/n > ^W-.^, 

As the function e~ 2 ^' x is concave down on x > /3/2, it follows that the left side of (|2ip is greater 
than 2(ie~P'^ n+l > /{n + l) 2 , the derivative of e -2 ^'* at i = n + 1, and it now suffices to show 
2 e 2/3/n e 2/3/(n+i) > ^ n _|_ i) 2 ^ 2 , p r this, simply note that the left side is always greater than 2, 
while (n + l) 2 /n 2 is less than 2 for n > 3. This gives the claim, hence w n < £e~ 2 ^' n for large n. 
In particular, we obtain v n <C £, hence u n < £X n /n for all n large. This completes the proof of 
Theorem 11.61 □ 
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Proof of Theorem \1.8[ We take all notation from Theorem 11.61 Note that our additional assump- 
tions imply that rk(n) = r^ is a positive integer for all k. We then have 

TV N 

/ 99 n _ V"^ a k v n-k ST~^ r k v n-k 

[ ' Vn+l ~ 2- ~>F^ + ^ x k +\n-k){n + l) ' 

fc=0 fc=o v j\ < J 

It follows that v n +i > min(i> n , . . . , i> n — jv) > for large n, by simply ignoring the righthand sum in 
([22]) . henceforth denoted R n +i- With Theorem 11.61 and the lower bound v n+ \ > min(u n , . . . , v u -n), 
this now implies that there exists an Aq such that < a < v n < £ for n > N%. By this, we 
have \R n \ < A/n 2 for some A > and n > N\. Fix an e > 0, and choose A^2 > N\ such that 
|i?n| < e/n 3 / 2 for n > N2. Now let M > A2 + A (M will be further specified later), and define the 
sequence {w n } n >M-N as follows: 

N 

w n +i = J2 \k+i ' Wj = Vj ioi j = M,...,M - N . 

k=0 

In other words, {w n } represents a linear approximation to {v n } (but with initial conditions that 
reflect the non-linear nature of the recursion for {v n } for small n). Next, we define 5 n := \w n — v n \, 
which satisfies the recursive inequality: 

N 

(23) 5 n+1 = \w n+ i - v n+ i\ < ^2 -rr^i&n~k + \R n +\\ < max(<5 n , . . . 5 n - N ) + — ^ 

fc=0 

n-M 

(24) < max(<5 M , . ■ ■ ,5 M -n) + £ ,.. ' W2 < C(3/2)e , 

kt^N ( M + k > 

by the fact that 5m -j = for j = 0,...,N. The manipulation of the maximum between (|23p 
and (J24l) is simply given by induction, using (|23l) . As the sequence {w n } satisfies a truly linear 
recurrence, we may express w n = Y^Pii 71 ) 00 ?' wherein p^ is a polynomial in n and O0i = Aj/A, 
corresponding to a scaling of the eigenvalues of x- We have 

K - %l < £(IM*)I + l«(v)l) • k ? " w "l ^ E (l«0»OI + l«G/)l) • (l«f I + k y D • 

By assumption, all Wj 7^ 1 have magnitude not equal to 1, and an application of Rouche's Theorem 
then implies that |wj| < 1 for oji ^ 1. Thus there exists an JV3 > such that x,y > A3 implies 
\wx — w y \ < e. Take M > A3. It then follows that 

\v n ~ Vm\ < \w n - v n \ + \v n - v m \ + \w m - v m \ < (2£(3/2) + 1)e , 

where we choose n,m > M. Thus {v n } represents a Cauchy sequence, and converges. Moreover, 
Theorem II .61 now gives linin^oo v n = £q = inf £, where the infimum is taken over all possible choices 
of £ in Theorem 11.61 Returning to the definition of u n , we now have u n ~ £o\ n /n, as desired. This 
completes the proof of Theorem 11.81 □ 

Suppose that the recursion in Theorems II. 61 and II .81 arises as the coefficient recursion to a Picard- 
Fuchs differential system. Then the condition a^ik — 1) + bk = is equivalent to the assumption 
that the Picard-Fuchs equation is of the form 

(25) A(t)F(t) + B'{t)F'(t) + B(t)F"{t) = A(t)F(t) + (B(t)F'{t))' = , 

with A(t),B(t) G Z[i] (with B(0) = 0). We note that the elliptic families associated to the 
congruence subgroups Ti(7), T(8; 4, 1, 2), and ri(10) correspond to differential systems of the form 
(I25D . For Ti(7) and ri(10), Theorem II. 81 applies directly. For T(8; 4, 1, 2), a similar result will hold 
under the condition that the associated sequence {v n } can be bounded away from (since Theorem 
11.61 does not apply). 
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Proof of Theorem M.KA Let vo = and v n = u n — a\ n /n. Since f(t) + olog(l — At) — > b as 
t — > (A -1 ) - , we obtain ^2v n t n — > b for t — > (A _1 )~. Our result in Theorem 11.61 implies that 
v n = 0(A n /n). Thus, Littlewood's extension of Tauber's Theorem [15j P-233-235] (and a change 
of variables) gives ^ v n \~ n = b, whence 

(26) ^u n A" n = alogx + 6 + 7 + o(l). 

n<x 

Let c(x) = u\ x \ A L^'J anc i define C{x) = ^c(t)dt. Since {u n X~ n } is eventually positive and 
monotonically decreasing, we have 

(s — r)xc(sx) < C(sx) — C(rx) < (s — r)xc(rx) 

for sufficiently large x and arbitrary real numbers r and s with < r < s. Since (|26p implies that 
C(sx) — C(rx) = a(logs — logr) + o(l) as x — > oo, it follows that 

logs -logr 

hmsupxc(sx) < a • 

x— >oo s r 

for any positive r < s. For s = 1 and r — >• 1~, we obtain lim sup^^^ xc(x) = a. Likewise, 

log s — log r 



liminf xc(rx) > a 



s — r 



for any positive r < s. With r = 1 and s —> 1 + , we obtain lim inL^oo xc(x) = a. Thus lim^^oo xc(x) 
exists and is equal to a, hence u n ~ a\ n /n. This completes the proof of Theorem 11.101 □ 

Theorem 11.81 do not allow us to determine £q exactly. If we are just working with a general 
holonomic recurrence, then the only option remaining might be to compute £o to a very high 
accuracy. However, if we are working with a holonomic recurrence which arises from a Picard- 
Fuchs differential equation, we can use Theorem 1 1 . 1 1 and explicitly find £q through modular forms. 
The referee kindly offered us the proof of Example 14.11 which applies Theorem 11.101 to show the 
asymptotic for the u n associated to Ti(7). The method to prove Examples 14.21 and 14.31 should be 
similar. 

Example 4.1. Consider the elliptic family associated to Ti(7). The equations parametrizing this 
family, the associated Picard-Fuchs equation, and associated recurrence were stated in (J3J) — dSJ) - 
With f(t) = Yln>o u nt n denoting the holomorphic solution of the Picard-Fuchs equation at t = 
and A ~ 6.295897 (the dominant root of x s — 5x 2 — 8x — 1), we show that f(t) + alog(l — At) —> b 
as t — > (A -1 ) - for some constants a and b which we give explicitly. 

We first identify the modular function t{r) and the modular form /(r) associated to the differ- 
ential equation ([5]). This can be done using the idea described on Page 1298 of [4]. We find 

E^rfE^T) E^rUrf 



E 3 (t)3 ' JW £!(t)3£ 2 (t)5 

in which 

oo 
E g (r) = q 7B (9/7)/2 "Q (1 _ /(n-l)+ g)(1 _ g 7n- 9) 



71=1 

„2 



denotes the generalized Dedekind eta function and B{x) = x — x + 1/6. 

To continue, we determine the analytic behaviors oft and / at various cusps of Ti(7), taken here 
as oo, 2/7, 3/7, 0, 1/2, and 1/3 (we need not consider the cusp class of 1/7). For this, we consider 
the actions of the matrices (J{), ( 7 I 3 ), ( 7 I 2 )> (? ~o)i (147); an d (217) on * anc ^ /• By the 
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transformation formulas for E g {r) given in |1Q|, Corollary 2], we have 
(27) t(r) < 2 "A - ^?E,{r) 



« -■>) E 6 (rf * + 

3 - 1\, _ E 3 (t) 2 E 6 (t) _ E^rfE^r) 

7 



(28) t(r) It - 2 J- ^ 9(T )3 =" ^ (T) 3 --1 + -- 

For the last three matrices, we note that E g (r) = i£ Sj o(7r), where E^o is defined in JT9], Page 673]. 
Now, let £ = e 2m ' . The transformation formulas given in |19l Theorem 1] imply 



E g (r) 



E 9 {r) 



E 9 (r) 



? "o 1 ) = E 9 ,o (=i) = e-(-V^/ 7 )^. 9 ( T ), 



7 


3 


14 


7 


7 


2 


21 


7 



2r + l, 

*i(2/3+V/7-3*/7) (_l)9 C -3 9 2 jE70)3g ( r ) ) 



£, 



9,0 



7r + 2 
3r + l 



=7 ri(-l/6+2 5 2 /7-2g/7) 



^7 5 ,2 9 (r) 



Then 



(29) 



t(r) 



t(r) 



(30) 



t(r) 



-1 

7 



7 3 
14 7 



7 2 
21 7 



3 27Ti/7 



e «(-l/6+2 9 2 /7-2 S /7) ( _ 1)SC -2 ff ^ o2ff(r) 

^o,-i(r) 2 ^o,- 2 (r) _ 



Eo,- 3 (r) Z 
sin 2 (67r/7)sin(27r/7) 
sin 3 (47r/7) 



+ 



^ (l- C -3)3 

0.15883360... + ••• , 



-6tti/7 



E , 3 (t) 2 E , 6 (t) 



E , 9 (tY 



sin 2 (4vr/7)sin(67r/7) 
sin 3 (2vr/7) 
7 E , 2 (t) 2 E 0A (t) 



+ 



-0.86293666 . . . + 



-iwi/ 



E 0fi (rY 



(31) 



sin 2 (27r/7)sin(4vr/7) 
sin 3 (67r/7) 



+ ••• 



-7.29589694... + 



In particular, we see that the singularity of the Picard-Fuchs equation at t « 0.15883360 corresponds 
to the cusp 0, i.e. A = l/i(0), in which A is the dominant root of x 3 — 5rc 3 — 8x — 1. Now 



V 



W 



r,(7r) 



whereby 



in which 



In addition, we have 



/-1\_ 7C 2 r^o,_ 3 (r)r ? (7r) 2 _ 
f \W) - ' £o,-i(r) 3 £o,-2(r)^ " ~ CZT{1 + '" h 



C 



7sin(47r/7) 



128sin 3 (67r/7)sin 5 (27r/7) 



n=0 



7rJ ' 
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which is valid for r is a neighborhood containing the upper imaginary axis. For convenience, let 
s(t) = t(— l/7r). Since (j~q) normalizes ri(7), it follows that s(t) is a modular function on Ti(7) 
and thus a rational function of £(r) of degree 1. Likewise, £ is a rational function of s(t) of degree 
1. To be specific, 

^ s(2/7) - 8 (3/7) s(oo) - s(r) 



s(oo)-s(3/7) s(2/7) -s(r)' 

because £(00) = 0, t(2/7) = 00, and £(3/7) = —1 (from our computation in ([27|) - ([28]) ). For r = ix, 
x > 0, our two expressions for /(— l/7r) give X]^=o u n s(ix) n = Cx(l + •••). Furthermore, we have 



-27ma:\(i„ 



E 3 (ixy ^x 

where d n = 0,2,1,-3 for n = 0,±1,±2,±3 (mod 7), respectively. Then x = — (log£(ix))/(27r) + 
o(l) as x — > 00, and it follows that 



v-^ C C ( s(ix) 

22u n s(ix) n = -— log£(ix) + o(l) = -—log II 



2-7T 2ir I s(oo, 

n=0 

- ^ log f s(oo)(s(2/7) - s(3/7)) \ 

2vr g V( S (oo) - s(3/7))(s(ix) - 8(2/7))) ^ l ' 
as a; — > 00. Therefore /(£) + alog(l — t/t(0)) — >■ 6 as £ — >• £(0)~, with 

a = — = 7sin(4^/7) ^ 0.3556270700876065 

2vr 256vr sin 3 (6vr/7) sin 5 (2vr/7) 

» = "- ^ f^ ^^(^l 377 ^ J - 0.7144010142820709, 
2tt S V(^oo)- S (3/7))( S (oo)- S (2/7)); 

wherein we have used the identities £(r)| ( - 1 °) = £(r), s(2/7) = £(-1/2) = £(1/2), and s(3/7) = 
£(—1/3) = £(1/3) as well as the numerical computations in (|29"l) - (i3~i~|) . At this point, an application of 
Theorem [TTTU] yields that u n ~ aX n /n, in which A ~ 6.295897, the dominant root of x 3 — 5x 2 — 8x — 1. 

Example 4.2. Next, consider the elliptic family associated to T(8;4, 1,2). From |16j we know that 
the defining equation is 

y 2 = x 3 - 2(8£ 4 - 16£ 3 + 16£ 2 - 8£ + l)x 2 + (2£ - 1) 4 (8£ 2 - 8£ + l)x. 

This gives a Picard-Fuchs equation of the form P 2 (t)f"+ Pi (t)f + Pof = where P 2 = £(£- 1)(2£- 
1)(8£ 2 - 8£ + 1), Pi = 80£ 4 - 160£ 3 + 102£ 2 - 22£ + 1, and P = 4(2£ - 1)(8£ 2 - 8£ + 1). Considering 
the associated recurrence for u n yields that 

2 (4 + 2^2)" 

u n ~ • 

7T n 

Example 4.3. Finally, we consider the elliptic family associated to ri(10) (an index 36 subgroup). 
Again from [16] we know that the defining equation is 

y = x(x + ax + b) 

with a = -(2£ 2 -2£+l)(4£ 4 -12£ 3 + 6£ 2 + 2£-l) and b = 16£ 5 (£-l) 5 (£ 2 -3£ + l). This gives a Picard- 
Fuchs equation of the form P 2 (t)f" + Pi (£)/' + P f = where P 2 = £(£ + 1) (2£ + 1) (£ 2 + 3£ + 1) (4£ 2 + 
2£-l), Pi = 56£ 6 +240£ 5 +320£ 4 + 156£ 3 + 12£ 2 -8£-l, and P = 2(36£ 5 + 128£ 4 + 136£ 3 +49£ 2 +2£-l). 
Considering the asymptotics of the associated u n gives that 

2 (1 + V5) n 

U n • 

^^5 + 2^/5 n 
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The constants which appear in front of the asymptotics for the recurrences associated to the 
congruence subgroups of index 24 and 36 are likely related to periods of the elliptic curves associated 
to the recurrence. Moreover, we conjecture that all the constants associated to such asymptotics 
will be an algebraic number divided by some half-integer power of n. 

5. Final Remarks and Future Work 

Systems of the form described in ([25]) have been studied by Beukers [3] and Zagier [20] due 
to their connection with Apery's proof of the irrationality of £(3) in [TJ; as such, their attention 
is restricted to the special case deg^4 = 1 and deg-B = 3. We briefly describe an extension of 
this system which arises in the case of certain index 24 congruence subgroups. Consider ()25p with 
instead degA = 3 and degi? = 5. Thus we have 

(32) (b 3 t 3 + b 2 t 2 + M + bo) F(t) + ((c 5 t 5 + c 4 t 4 + c 3 t 3 + c 2 t 2 - t)F'(t))' = . 

If 63 = C5 = 0, then the associated holonomic recurrence will have order 3. These we choose to 
ignore, as none of the recursions in the index 24 subgroups are of order less than 4. Under the 
parameter change t h-> —t, our coefficient series undergoes u n h-> (— l) n u n , hence we may suppose 
C2 > without loss of generality, and concentrate a brute force search on the domain of about 28.8 
billion 8-tuples given by 

\h\ < 10, l&al < 25, |&i| < 15, |&o| < 5, |c 5 | < 5, |c 4 | < 10, |c 3 | < 15, and < c 2 < 10. 

Searching through seems to give evidence to a variant of one of Zagier's conjectures in [2D]; that is, 

Conjecture. Any integral solution to the differential equation in (]32p corresponds to the solution 
of a Picard-Fuchs equation about a singular fiber. 

While in this paper we concentrated on index 24 subgroups, specifically ri(7), the authors expect 
that the same analysis can be applied to ri(10) and other subgroups of index 36, 48, and 60 though 
such analysis might be more technical. 
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